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Abstract 

We derive a new class of exact time dependent solutions in a warped six dimensional supergravity 
model. Under the assumptions we make for the form of the underlying moduli fields, we show that 
the only consistent time dependent solutions lead to all six dimensions evolving in time, implying 
the eventual decompactification or collapse of the extra dimensions. We also show how the dynamics 
affects the quantization of the deficit angle. 
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I. INTRODUCTION 



Six dimensional supergravity models have several interesting properties. Salam and Sez- 
gin obtained static solutions in which the six dimensional gauged supergravity compacti- 
fies on a product spacetime of four dimensional Minkowski and a two dimensional sphere, 
M.a x S 2 Remarkably this supergravity model admits the supersymmetric Minkowski 
vacuum while many other supergravity models do not. The modern interpretation of this 
property is that the solution of this theory is compatible with the introduction of branes 
into the spacetime. As with any massive defect, this then leads of course to the appearance 
of a deficit angle in the two internal spatial dimensions, as a gravitational response to the 
tensions of the branes 2|, y) . The resulting geometry looks like a rugby ball solution where 
the branes are located at the north and south pole of the ball. Gibbons, Guven and Pope 
(GGP) [4] showed that the Salam-Sezgin vacuum is in fact the unique one with a four di- 
mensional maximal symmetry and general static solutions with an axisymmetric internal 
space. 

The observation in Ref. {4] that the four dimensional spacetime is always Minkowski even 
in the presence of branes with tensions forms the basis of the interesting supersymmetric 
large extra dimension (SLED) scenario, a recent approach to solving the cosmological con- 
stant and dark energy problems j^]. If only for this reason, such is the prize at stake, it 
makes this six dimensional supergravity interesting from a cosmological point of view, al- 
though we note that cosmology in six dimensional supergravity has previously been studied 
in the context of Kaluza-Klein cosmology 

HQ. 

One of the neatest aspects of the SLED 
model is that a 3-brane with any tension in six dimensional spacetime induces only the cor- 
responding deficit angle and maintains a vanishing four dimensional cosmological constant, 
at least at the classical level. This feature is often referred to as a "self-tuning mechanism" 
of the effective four dimensional cosmological constant and would be expected to be part of 
the solution to the cosmological constant problem (although we still have to account for the 
affect of quantum corrections). 

Although the SLED scenario has enjoyed a number of successes, open questions still 
remain. We are particularly interested in establishing whether the self-tuning mechanism 
really works in a time dependent evolving Universe (another attempt, see |7|). Previous 
authors have argued that the self-tuning of the four dimensional cosmological constant does 
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not work at least in non-supersymmetric six dimensional Einstein Maxwell theories 







The work we present here has a number of overlaps with that of Tolley et al 
12j |. (and more recently with that of Kobayashi and Minamitsuji [3]) who have obtained 



a series of solutions to the dynamical system based on an elegant scaling argument. We 
believe that the explicit expressions presented here of exact time-dependent solutions to the 
six dimensional supergravity model are given for the first time. We begin by deriving the 
underlying field equations in section [Til This is followed in section [TTT1 with a demonstration 
that it is impossible to have static internal spaces with a corresponding expanding external 
three dimensional space. We extend the analysis to a fully time dependent case in section 
[TV] and find a new class of exact solutions showing the nature of the instability and resulting 
evolution of the compact dimensions. Finally we conclude in section [Vj 



II. THE BASIC FIELD EQUATIONS 

Concentrating on the bosonic field contents of this model, we have the metric gMN, dilaton 
4>, a U(l) gauge field A M with field strength F MN and an antisymmetric tensor field B MN 
whose corresponding field strength is expressed as 

Gmnp = QmBnp + FmnAp + cyclic permutations. (1) 

The lagrangian density for the bosonic sector is given by 

£sugra = \R- \d M <pd M <P - e ^-G MNP G MNP - ^F MN F MN - 2<M (2) 

where g is the U(l) gauge coupling 1 . Here M,N run over all the spacetime indices and 
we work on the two - sphere of radius r with the six dimensional (reduced) Planck scale 
M 6 = l. 

The field equations are 

□0 + —^—G mnp Gmnp + ~7~ F MN Fm n - 2g 2 e* = 0, (3) 
o 4 

D M (e- 2 *G MNP ) = 0, (4) 
D M ( e ^F MN ) + e~^G MNP F MP = 0, (5) 

1 Note that our definition of cj> and g are different from those in GGP. They are related through — 2(j) OUYS — 
0ggp and gl UIS = 2g^ GP . 
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—Run + dM<fidjy4> H — — ^Gm PQ Gn P q — -G OPQ GopQgMN S j 

+e~* (f m p F np - ^F p ^F PQ g MN ^ + g 2 e*g MN = 0, (6) 

and following the usual ansatz adopted for simplicity, from now on, we consider the case of 
a vanishing three form field strength 



qMNP = Q (?) 



Then, above field equations become 



□0 + 6 —F MN F MN - 2 5 V = 0, (8) 
D M ( e -*F MN ) = 0, (9) 
-Run + d M (j)d N (j) + e"* (f m p F np - ±F PQ F PQ g MN ^ + g 2 e^g MN = 0. (10) 

The metric ansatz we adopt is 

ds 2 = U{x m ,t) 2 ds A + r{t) 2 ds 2 2 , 

ds\ = -dt 2 + 5 ij dx l dx j , ds\ = j mn (x m )dx m dx n , (11) 

where i, j run over the usual three spatial indices, m, n run over the extra two spatial indices 
and r jmn{x m ) is an arbitrary two-dimensional metric. With this metric, the two form field 
strength takes the form of 

F = F = 

Fmn F(tjX )6 mn , (12) 

with e mn being the anti-symmetric tensor. Here, the use of Greek indices denote the four 
spacetime coordinates (i.e. n = (0, i)). With the metric ansatz Eqn. (fTTj) . we write each 
component of the Einstein equations, (0 — 0), (0 — m), (i — j), (m — n) respectively as 

2" 



.0 



d r\ d rd U ( d r 



r j r U \ r 



-y 2l mn d m Ud n U - ^D m (^dnU 2 ) + e -^F PQ F PQ U 2 - g 2 e <t >U 2 = 0, (13) 

d ™U\ ( d U\ d m U d r 

+ 4 — — - 4— + a o 09 m = 0, (14) 



U J V U J „ C/ r 



U J V U J U r 



+2 drnUdnU^ mn + J-D^mnQ^ _ tlpPQp jfl + gW = q (15) 

U U r z 2r z 8 



d r 2 \ ^d U d r r 2 



2U 2 J U r U 2 



u i - u u ' 



(16) 



+d m <f>d n <f> + e~* ^F m p F Pn - ^F PQ F PQ g mn ^ + g 2 e <t> g mn = 0. 
A. The static Gibbons, Guven and Pope solution 

Before discussing the time dependent solutions of this system, we recall the derivation of 
the static solution originally obtained by Gibbons, Guven and Pope (GGP) [4]. When we 
take r(£) = 1, U(x m ,t) = W(x m ) and (j) = <p(x m ), the Einstein equations and the equation 
of motion for the dilaton <fi reduce to 

1 D m { 1 mn d n W 4 ) + e -^F P( ^F PQ -g 2 e^ = Q 1 (17) 



4W 4 



\ W J WW 

-d m (j>d n (f> - e 4, (^F m p F nP - ^F PQ F PQlm }j - g 2 e^ mn = 0, 



(18) 



and 



^D m (V^ 4 ^) + 6 -^F MN F MN - g 2 e* = 0, (19) 
leading to the GGP solution: 

0=-^lnJ¥ 4 . (20) 



The equation of motion (J9J) for F mn leads to the solution to Eqn. ([12]) 

F = -|e*iy- 4 = ~W~ 6 , (21) 

where we have used Eq. (TSUI) . Note that q can be interpreted as a magnetic charge. 

It will prove useful to recall the explicit solution for <fi and j mn as presented by GGP 

0,Q 

( e -<f> Q-<t>r 2 \ 

e 2 * = (23) 
h 
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with 



/o = l + 3. A = l + 3. ( 24 ) 
'o ' 1 

r 2 = A> ^ = 4- (25) 

In the following sections we begin to explore the dynamical equations by allowing the 
scale factors and the fields to become time dependent. 

III. ONLY ONE EVOLVING SPACE IS NOT A SOLUTION 
A. Time dependent U 

Ideally what we want to obtain is a solution which describes the expansion of our three 
space dimensions with a static extra dimensional space. As a first step towards obtaining 
it, we make the ansatz of a static internal space r = 1, and a static dilaton (p = <j)(x m ). This 
combination makes sense in that the relation r 2 oc has previously been obtained in 
hence the dilaton would be static if r is static. In fact, in the following subsection, we will 
show that allowing for a time dependence of cf> does not improve the possibility of obtaining 
static r solutions. The field equations are reduced to 

§S (ir) q " 4^4 D m (l mn dnU*) + ^-F PQ Fpq - ,V = 0, (26) 

^ =0, (27) 



TP &) o " (^) 2 ^ + - X^** + ^ = °< <»> 



^DA d -f-)+^ d -f-R m ' 



+d m (f>d n (f> + e~* (F m p F Pn - }-F PQ Fp Qlm }j + g 2 e^ mn = 0, (29) 
j- 4 D m (r^d n f\ + 6 -^F MN F MN - «?V = 0. (30) 



Although the equations look intractible, we can make progress by noting that because 

d m U\ d d m U d Ud m U ( d U' 



u n U U 2 V u 



(31) 
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then from Eq. ([2 

'd m U\ ( d U 



U J, U 



0. (32) 



The general solution of U is therefore given by \nU = lna(t) + lnW(x m ) where a(t) and 
W(x m ) are integration functions. Thus, we find that U has to take the separable form of 
U = a(t)W(x m ). The field equations f[2"B"l) and f[2"8l can then be reduced to 

C(x m ) - J^D m (j mn d n W 4 ) + e -^F PC ^F PQ - <?V = 0, (33) 
3 / doa\ 1 / doa\ 2 / <9oa x 2 



.0 

which after some algebra leads to 



, > i / / m i + 772 — =C(x m ), (34) 



a(t) = (35) 
C(x m ) = 2T J (36) 

with a and t being integration constants. 

However, from Eqs. ( J26l) and ( J30l) . we also know that 

An ( 7 mn f/ 4 ^(ln f/ 4 + 20)) - 4f/ 4 C(x m ) = 0, (37) 

which is conflict, for any non- vanishing U on the internal space, with the fact that the extra 
two dimensional space is compact. In other words, if we integrate both sides of Eq. (I3T|) 
over the compact extra space, we see that the first term on the left hand side vanishes as it 
is a total derivative while the second term does not. Hence, there is an inconsistency and so 
we conclude there is no static solution for the compact space with this ansatz for U. 

There is a caveat to this argument. We have implicitly assumed that the extra space is 
smooth. However, if we allow the extra space to be singular then it is possible that, de-Sitter 

n 

type solutions may be obtained |16J. 

B. Time dependent U and 0, but static r. 

We now allow for a time- dependent dilaton. The field equations are 

3 (j^j q + do W - ^D m (r n d n U*) + e -^F p QF PQ U 2 - g^U 2 = 0, (38) 
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+ 4 f^l + d <j>d m <j> = 0, (30) 



u J >0 \u / MI 

2 



(ir) o + 2 (ir) ~ W- Dm ^ mnQ ^) + 6 ^f pq FpqU 2 - *W» = 0, (40) 



.0 

+d m <R</> + e"* ^F m p F P „ - ^F PQ F PQ g m }j + 9 2 e^„ = 0, (41) 
~j- 4 do (u%fj + ijD m (V 4 7™^4) + ^F»»F MN - 9 2 e* = 0. (42) 

If we again assume the form U = a(t)W(x m ), then the (0 — m) component of the Einstein 
equation implies <f> = <p(t). In addition, from the equation of motion, using the solution for 
the flux F oc U(t, x m )~ 4 , we find that W = 1 because each term has a different dependence 
on W and the cj) equation of motion can not be satisfied if U depends on x m . Hence we must 
have U = a(t) and moreover since the spacetime no longer has a non-trivial warp factor 
W(x m ), it can not be warped. Given the above result, the equations of motion now can be 
written as: 

3 (^j + + ^-^-F PQ F P qOl 2 - gW = 0, (43) 

~h 9 ° 2 ) + e -Y FMNFMNa2 ~ 9 W = °< ( 45 ) 

Notice that the term g mn R m ' mm 'n could in principle be a function of x m , but in this case it 
is not allowed by Eq. (146p as the left hand side depends only on t. This fact means that the 
compact extra space must be a constant curvature two dimensional sphere. Here there is no 
way to introduce branes which induce a deficit angle and deform a sphere with a constant 
curvature into a rugby ball shape. Therefore, we can see that this ansatz, namely varying 
U and (f> with static r can not lead to satisfactory solutions. 
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C. Time dependent r and (f>, but static U. 



Finally, let us try to obtain a solution of the static three space with a dynamical extra 
dimension. If we take r = r(t) and (ft = <p(t, x m ), but U = W(x m ), then Einstein's equations 
and the equation of motion for <f) are 

+^-^^ + ^-W^0, ,47, 

-4%^ + S„0S„0 = 0, (48) 
W r 

+a m <R0 + e~* ^F m p F Pn - ^F PQ F PQ g mn ^ + g 2 e*g mn = 0, (50) 

' ;d ( r*d ^) + J^D m ( W^dnt) + ^F MN F MN - = 0. (51) 



W 2 r 2 »y - 2 J W 4 r 2 V 2/8 

Now provided that x m ) can be decomposed as 0(t, s m ) = </>(£) + (f}(x m ) and F depends 
only on x m , then Eqs. ()47l) and (1491) can be reduced to 

r(t)V (t) = 1, (52) 



,0 



71 (^ mn fj W 4 ] P -<t>{x m ) 



(54) 



Unfortunately, the solution of these equations are not compatible with 

d r 2 = 0, (55) 

which can be obtained from Eq. (l5Tj) using Eq. fl52|) . Thus, we once again see that there is 
no consistent solution with this ansatz. 

IV. TIME DEPENDENT SOLUTIONS WITH DYMANICAL r, U AND <f> 

Having tried unsuccessfully to obtain static solutions for r and U, we now look for dy- 
namical solutions where all the key fields r, U and <f> are time dependent. We again make 
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a series of ansatz, in this case (p(t,x m ) = <fi(t) + 0(x m ), and assume the separable form of 
U = a(t)W(x m ). The Einstein equations and the equation of motion for are: 



3 / d a 
U 2 \~ 



+ 



U 2 



d r 



dordoa | / d r 
r a \ r 



W u o<, 

u 2 



1 D m ( 7 m "cW 4 ) + ^-F P *F PQ - <?e* = 0, 



4r 2 W 4 
dmWdpr 

W r 

1 / d a 

U 2 \ a 



0. 



(56) 
(57) 



J2_/<9oa\ 2 d ad r 
U 2 \ a U 2 a r 



4r 2 W 4 



' D m (j mn d n W*) + ^F P ®F PQ - g*e* = 0, 



(58) 



W 7 



w w 



( d r 2 \ ^ dpadpr r 



\2W*J fl 

+ e~+ (F m p F Pn - ^F PQ F PQ g m ^j + g 2 e^g mn = 0, 



a r U 2 



7n 



(59) 



-J^rA (W^) + ^D m (wS™d n ^ + ( —^-F mn Fmn - *V = 0. (60) 

Jnder the additional ansatz that e^*V 2 = 1 (motivated by the observation that r 2 oc 
3|), which is equivalent to 



d (, 



d r 



we see that Eq. ( 1571) . leads to 



0„ 



-2- 



w ' 



(61) 



(62) 



as in the GGP solution. If we further assume that the field strength F is static and only 
depends on x m , then the field equations Eqs. (15B1) . (158"]) and (I6"0l coupled with Eq. (I6"2"j) can 
be rewritten as the following differential equation which depends only on x m , 



(63) 



where, C(x m ) is given by 



r 

U 2 



d a 



d r 



/" ■'"■> .2 ( iir ) a 

r a \ r 



d a 



,' dna\ dnadnr 
+ 2 ( — + 2— — — 
a I a r 



each equality in Eq. (164|) arising from Eqs. fl56|) . (1581 and (160]) . respectively. 



(64) 
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A. Power law solutions for r and a 



Eqns. (1631) and (I64|) still look very difficult to solve directly from first principles, and so 
instead we will try to obtain solutions by assuming the form of r and a, and looking for 
self-consistency in the solutions. As a first attempt we assume power law behaviour for 
them, namely: 

a<xt n , roct nr . (65) 
The three equalities in Eq. (1M1) coupled with Eq. flBTj) now become 

-j2(n r — n— 1) 

C{x m ) = — — — (-3n - 2n r - 2n r n + 6n 2 r ) 

-j2(n r —n—l) 

= 777^ n r (2n + 2n r -l) 

W 2 

= — 2 n{2n + 2n r -l). (66) 

There are two possible ways in which C(x m ) can be a function of only x m , as required 
by Eq. (163]) . The first is if the time dependent prefactor vanishes which corresponds to 
n r — n — 1 = 0. The second way is if the right hand side of each of the terms vanish, which 
corresponds to the brackets vanishing in Eq. (1661) . The former is precisely the structure 



found by Tolley et al 



12] based on a scaling argument for the scale factors. However, the 



metric ansatz of g^ v in 12j is slightly different to ours ([ITT) . In particular it follows that the 
condition n r — n — 1 = is not a solution in our case, because it can not satisfy all three 
equalities in Eqs. fl66l) . In fact we determine the values of n and n s in Eqs. f l66|) by equating 
the coefficients : 

— 3n — 2n r — 2n r n + 6n 2 r = n r {2n + 2n r — 1) = n{2n + 2n r — 1). (67) 
This has the non trivial solution 

2±y/3 _^V3 
n= — , n r = T—- (68) 

This in turn gives C(x m ) = which is consistent with the above discussions and is compatible 
with Eq. ( !59l) too, because the solution satisfies 

d r*\ +A doad«rr* =Q (gg) 



2a 2 



n a r a? 
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Notice that in this case, we obtain identical solutions for F, <j) and W as found in the GGP 
solution. This is as expected, since the x m dependent part of the field equations are identical 
to that of the GGP solution. In this sense we have obtained the time dependent version of 
the GGP solution. 

From the metric ansatz Eqn. fill I) it follows that the time t is actually the conformal time 
in the usual sense. The "cosmic time" r can therefore be defined as dr oc t n dt, from which 
we obtain 

ds 2 = W{x m ) 2 [-dr 2 + a(r)%cfaW] + r(r) 2 ds 2 2 , 

a(r) oc r n/(n+1) , r(r) oc W (n+1) , (70) 

in terms of the cosmic time. 



B. Exponential solutions for r and a 

The next obvious step is to assume an exponential form 

a(t) = e ht , r(t) = e M , (71) 
where h and h r are constants. In this case, Eq. (1641 with Eq. ([ST]) leads to 

„2{h T -h)t 

C(x m ) = w2 2h r (-h + 3h r ) 



e 2(h r -h)t 

w 2 

e 2(h r -h)t 



2h r (h + h r ) 

2h(h + h r ). (72) 



W 2 

which now has a non trivial solution h = h r . Then, C(x m ) is given by 



Thus, we obtain the equations of motion of the x m dependent part of the fields to be 

^pq ~ 9 2 e^ m) = 0, (74) 



Ah 2 T) ('y mn r) W A \ P -<t>( x,n ) 

\ n u m [pi a n w ) + e F p Q _ 2 e ^) = 



W(x m ) 2 AW 4 
Ah 2 



W{x™) 2lmn ~ dm<t>dn<t) ~ 6 ~ HXm) ( FmF>Fpn ~ I fPQf pq1™?J - g 2 e^ xm) lm n = 0.(75) 

Something significant can now be seen. Recall that we have equation (1621) . relating (f) and 
W(x m ). Given the solution we have just obtained, we see that in Eqn. (1751) . by introducing 
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g 2 = g 2 — Ah 2 , then a new solution to the system is obtained which looks identical to the 
original x m part of the GGP solution but with our redefined gauge coupling g 2 replacing 
the original g 2 coupling. Obviously, the /i — *> limit corresponds to the original static GGP 
solution. Hence, we have obtained the explicit expression of the solution including the x m 
dependent part. However, notice that since h is just a constant it could in principle take 
any value. In particular, for Ah 2 > g 2 , corresponding to a negative g 2 , we find that the x m 
dependent part of the solution has only differs slightly from that obtained in GGP. We show 
this and give the actual solution for the case of vanishing and negative g 2 in Appendix B. 
The line element of this solution with such a nonvanishing h is rewritten as 

d s 2 = W{x m ) 2 [-dr 2 + (hT) 2 S ij dx i dx i ] + (hr) 2 ds 2 2 , (76) 



121 ]. however, 



in terms of the cosmic time. This solution is the same as that found in Ref. 
the x m dependence of the solution was not explicitely solved for there. Here, we have shown 
that it is same as that of the GGP solution. 



V. CONCLUSIONS 

We have derived a new class of exact time dependent solutions in a six dimensional 
gauged supergravity compactified on a two dimensional axisymmetric space. Under the 
assumption of a separable form of U we showed that there is no solution expressing the 
either an expanding four dimensional universe with a static internal space or visa versa. 
Exact solutions we obtained involved all the dimensions either expanding or contracting 
which means the eventual decompactification or collapse of the extra dimension, indicating 
an instability of Salam-Sezgin, (Minkowski) 4 x S 2 , spacetime for the case with the absence 
of the maximal symmetry in the four dimensional spacetime. 

In the above analysis, we did not include into the action brane terms such as 

Sbrane = ^ / ^ T i = J2 f ^ ^ ^ ~ > ( 77 ) 

i i 

where Ti is the tension of the 'i-th' brane and x™ denotes the position of the brane in the 
internal space. However, we can easily introduce such brane terms, their affect being to 
induce the deficit angle in the internal space. The topological condition for the gauge field 
A M is the same as we previously obtained for the static solutions, because the solution of 
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;he gauge field strength F mn is unchanged in the presence of the branes. As is discussed in 
4j, for the case of r 2 ^ r\ in Eq. fl24|) . while one pole can be smooth, the other has a deficit 
angle 

£ = 1-4- (™) 

27T Tq 

Combining Eq. (1251) and the topological condition, the Dirac quantization condition, for the 



gauge field Am becomes 



13, 



*t = N, 

q 



(79) 



leading to the quantized deficit angle 



2tt 



1 - N 2 



(80) 



which was previously obtained for the static solution with N being an integer J4]. However, 
for the new solution in Section IV B, Eq. ( 1751) is rewritten as 



i- = i-3£ 



2tt 



for a positive g , 



for a vanishing g 2 and, 



A = 1 

2tt 



8(~^ 2 ) 
q 2 ' 



for a negative <jr . Hence the deficit angle is given as 



2tt 



1 - N 2 



1-4- 



^ 2 



(81) 



(82) 



(83) 



(84) 



This implies that the interval of the quantized deficit angle becomes narrow for 2h ~ g in 
the time-dependent solution. It would be interesting to investigate the consequence of this 
new deficit angle. 
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APPENDIX A: CONVENTIONS AND GEOMETRICAL QUANTITIES 



Here, we note our conventions and several geometrical quantities. 
Christoffel symbols 

r 7VP = ^9 MQ (9qn,p + 9qp,n - 9np,q)- 



(Al) 



Riemann tensor 



-ft NOP — J- Arp,0 ~~ 1 JVO.P 1 QO 1 ATP — 1 QP L NO- 



(A2) 



With this definition, the sign in front of the Einstein term in the action is a plus. 
Metric 

ds 2 = U 2 {t,x m ){-dt 2 + 5 l3 dx l dx j ) + r(t) 2 1 {x m ) mn dx m dx n , 



(A3) 



where run over the usual three-spatial dimensions and m,n,..., run over the extra 

spatial dimensions. 
Christoffel symbols 



pO _ 1 d U 2 
1 00 — 2 7J 2- 

pO _ Idyll 2 s 
1 ij ~ 2 U 2 °V 

-pi _ 1 dpU 2 Xi 
1 0j - 2 U 2 3 

rra 1 d n U 2 n ,mn 

00 2 r 2 ' 



TO _ 1 9mC/ 2 

1 Ora — 2 I/ 2 

pO _ \ihr± 

1 mn ~~ 2 t/ 2 ^ mn 

pi _ 1 d m U 2 xi 
1 jm 2 C/ 2 J 



On 2 r 



P m — 1 d n U 2 mn x rm _ l^mqf^ i ~ _ \ 

1 ij 2 r 2 ' *J np 2 ' V '9 n >P ~ '"Pi?/ 

others = 



Ricci tensors 



00 



d U 
U 



-7™<9 m [/d n [/(<^ - 1) 



.0 



+ 



r /,o r V r 



ft 



0m. 



m CA _ fdoU 



R 



U 

d U 
U 



(S\ + 1) 



u 



+ 



d m ll d r 



D m h mn d n U 2 ) 



,m 



U r 



2r 2 

(5\ + S\-l), 



+ 



,o 



d u\ k d m ud n uu 2 k 

{d k -l) -7 (d k - I) 



U 



U U r 2 



d Q Ud Q r D m (l mn d n U 2 ) 

H 77 m 



U r 



2r 2 



(A4) 
(A5) 



(A6) 
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Rr. 



(S\ + l)D 
dor 2 



d m U 



2U 2 



U 

d U d r r 2 



(5\ + 1) 



d m U d n U 



,o 



2 „2 



others = 0. 



(A7) 
(A8) 



APPENDIX B: SOLUTIONS FOR VANISHING AND NEGATIVE g 2 

In this Appendix, we note the solutions for the case of vanishing or negative g 2 . Following 
GGP, by introducing the variables 

x = -(f) + In A, 

y = ^0 + 4lnW + In A, (Bl) 
z = -(f) - 2 In W, 

we obtain 

<v 1 



-in) +A 9 2 e 2y = \l (B2) 

ri 

and similar equations for x and z, both of which are decoupled from y Here, A2 is a 
constant for the first integral and rj is a coordinate in the coordinate system 

ds 2 2 = W s A 2 drf + A 2 dtp 2 . (B3) 



The solution for a positive g 2 is presented in Ref. J4] as 



1. / AI 



|/= -lncosh(A 2 (r/-?7 2 )) + -ln^^yJ . (B4) 

However, as one can see, in Sec. IV B, the effective g 2 , namely g 2 , can be positive or 
negative in some time dependent solutions. This then means that the solution for y is 
replaced with 

y = X 2 (v-V2), (B5) 

for the case of vanishing g 2 and 



\2 



y = -msinh(A 2 (77-77 2 )) + ~ln( ) , (B6) 

for a negative g 2 . 
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We then obtain 



AW 



and 



1 cosh(Ai(?7-?7i)) / Aj 

(cosh 3 (Ai (77-771 )) cosh(A 2 (r?-»?2 ))) 1/4 cosh(A 2 (77-772)) \4g^ 

1 cosh(Ai (77-771)) / Aj 

(cosh 3 (Ai(?7-7?i))smh(A2(??-??2))) 1/4 sinh(A 2 (77-772)) v -4<? 2 



1 



cosh(Ai(7?— 



(cosh 3 (Ai(?y-7yi))e- A 2('7-'72))i/4 e A 2 (77->ja) 



2Af 



-1 



for positive g 



^2 ) for negative g 



for vanishing g 



9 2 

7,2 



(B7) 



A 



(cosh 3 (Ai (97-171 )) cosh(A 2 fa-Tfe))) 1 / 4 V 4 ^ 



1/2 



(cosh 3 (Aj (77-771 )) S inh(A 2 (77-772 ))) 1/4 



Al V'' ( I 2 



v 1/2 

V 2At 

-3/2 



-3/2 



-3/2 



7,2 



for positive g 
for negative <? 2 
for vanishing ^ 2 



:bs) 



k (cosh 3 (Ai(r ? -r ? i))e- A 2('7-'72))i/4 \^2A^ 

Here, Ai is a constant for the first integral with repect with x, and we used the solution of 
x given in 4j. Setting Ai = A 2 = 1 and introducing a new coordinate dr = AW 4 drj, allows 
us to derive the deficit angle given in Eqs. (IHTj) . (|82j) and (1551) . 
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